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K — >• 7T7T amplitudes from lattice QCD with a light charm quark 
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We compute the leading-order low-energy constants of the AS = 1 effective weak Hamiltonian in 
the quenched approximation of QCD with up, down, strange, and charm quarks degenerate and light. 
They are extracted by comparing the predictions of finite volume chiral perturbation theory with 
lattice QCD computations of suitable correlation functions carried out with quark masses ranging 
from a few MeV up to half of the physical strange mass. We observe a AI — 1/2 enhancement in 
this corner of the parameter space of the theory. Although matching with the experimental result 
is not observed for the A7 = 1/2 amplitude, our computation suggests large QCD contributions to 
the physical A7 =1/2 rule in the GIM limit, and represents the first step to quantify the role of 
the charm quark-mass in K — > nn amplitudes. The use of fermions with an exact chiral symmetry 
is an essential ingredient in our computation. 

PACS numbers: 11.30.Rd, 12.38.Gc 13.25.Es 



INTRODUCTION 



The decay of a neutral kaon into a pair of pions in a 
state with isospin / is described by the transition ampli- 
tudes 



iA ie t5 ' = ((vr^/liJwIif ), I = 0, 2, 



(1) 



where H w is the AS — 1 effective weak Hamiltonian and 
5 1 is the 7T7r-scattering phase shift. The well-known ex- 
perimental fact 



\A /A 2 



22 



(2) 



is often called the AI =1/2 rule. Many decades after its 
experimental discovery, it is embarrassing that the origin 
of this enhancement is still not known. In the Standard 
Model (SM) a reliable perturbative computation of short- 
distance Quantum Chromodynamics (QCD) corrections 
[1-4], together with a naive order-of-magnitude estimate 
of long-distance contributions, would suggest comparable 
values for \Ao\ and \A 2 \ [1, 2). The bulk of the enhance- 
ment is thus expected to come from non-perturbative 
QCD contributions, which makes the AI — 1/2 rule one 
of the rare cases where an interplay between strong and 
electro-weak interactions gives an opportunity for a re- 
fined test of non-perturbative strong dynamics. 

Lattice QCD is the only known technique that allows 
us to attack the problem from first principles and possibly 
to reveal the origin of the enhancement [5, 6]. It would 
be interesting to understand whether it is the result of 
an accumulation of several effects, each giving a moder- 
ate contribution, or if it is driven by a dominant mecha- 
nism. Recently we proposed a theoretically well defined 
strategy to disentangle non-perturbative QCD contribu- 
tions from the various sources [7] , and in particular to re- 
veal the role of the charm quark and its associated mass 



scale (whose relevance in this problem was pointed out 
in Refs. [8, 9]). The main idea is to compute the leading- 
order low-energy constants (LECs) of the CP-conserving 
AS = 1 weak Hamiltonian of the chiral low-energy ef- 
fective theory as a function of the charm quark mass. 
They can be extracted by comparing finite-volume chi- 
ral perturbation theory (ChPT) predictions for suitable 
two- and three-point correlation functions with the anal- 
ogous ones computed in lattice QCD at small light-quark 
masses and momenta. The suggestion of using ChPT in 
connection with kaon amplitudes was pointed out long 
ago [10, 11]. It is only now that these ideas can be formu- 
lated and integrated in a well defined strategy [7] , follow- 
ing significant conceptual advances in the discretization 
of fermions on the lattice as well as enormous gains in 
computer power. The main theoretical advance is the dis- 
covery of Ginsparg- Wilson (GW) regularizations [12-14], 
which preserve an exact chiral symmetry on the lattice at 
finite lattice spacings [15]. By using these fermions the 
problem of ultraviolet power divergences in the effective 
Hamiltonian H w [16] is avoided in the case of an active 
charm [17], and quark masses as low as a few MeV can 
be simulated. Eventually the full K — > 7T7t amplitudes 
can be computed using finite- volume techniques [18, 19]. 

The aim of this letter is to report on a computation of 
the LECs of the CP-conserving AS = 1 weak Hamilto- 
nian with up, down, strange, and charm quarks degener- 
ate and chiral (GIM limit), i.e. the implementation of the 
first step of the strategy proposed in Ref . [7] . We perform 
the first quenched lattice QCD computation of the rele- 
vant three-point functions with quark masses as light as 
a few MeV, which turns out to be essential for a robust 
extrapolation to the chiral limit. Our results reveal a 
clear hierarchy between the low-energy constants, which 
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in turn implies the presence of a AI = 1/2 enhancement 
in this corner of the parameter space of (quenched) QCD. 

Since we are looking for an order-of-magnitude effect, 
and since simulations with dynamical fcrmions are very 
expensive, it is appropriate for us to first perform the 
computation in quenched QCD. The latter is not a sys- 
tematic approximation of the full theory 1 . However, 
when quenched results can be compared with experimen- 
tal measurements, discrepancies of 0(10%) are found 
in most cases [21]. In the past there were several at- 
tempts to attack the problem by using quenched lattice 
QCD [22-27]. In particular, in Refs. [25, 26], a fermion 
action with an approximate chiral symmetry was used 
and, despite the fact that the charm was integrated out 
and therefore an ultraviolet power-divergent subtraction 
was needed, the authors observed a good statistical signal 
for the subtracted matrix elements in a range of quark 
masses of about half the physical strange quark-mass. 
Several computations of Aj which use models to quantify 
QCD non-perturbative contributions in these amplitudes 
can also be found in the literature (see Refs. [28, 29] and 
references therein). 

THE AS = 1 EFFECTIVE HAMILTONIAN 

In the SU(4) degenerate case and with GW fermions, 
the CP-even AS = 1 effective Hamiltonian is [1, 2, 7] 

where 

Qf = Z±{(s^P-u)(u^P-d) 

± (s^P.d)(u^P.u) -[«-<;]}, (4) 

and any further unexplained notation in the paper can 
be found in Ref. [7]. We are interested in the ratios of 
correlation functions 

d±/ \ C^(x ,y ) 

where 

C(x ) = ^<[Jo(x)] Q/3 [Jo(0)W, (6) 

X 

Ct(x Q ,y Q ) = J^([Jo(x)}du[Qt(0)}[Jo(y)}us) , (7) 

[J/j]ap — Zj('4>aJnP~ipf3), and Zj is the renormalization 
constant of the local left-handed current. 



On the other hand the ambiguity in the definition of the LECs 
pointed out by the Golterman and Pallantc [20] is not present in 
the GIM limit. 



In the chiral effective theory the corresponding effec- 
tive Hamiltonian reads 

= ^§rV: s V ud {gtQt+9iQl} , (8) 
where, at leading order in momentum expansion, 

Qf = ^-{(ud^) us (uo^) du 

± (Ud^) ds (Ud^) uu -[ u - c ]}. (9) 

The complete expressions at the next-to-leading order 
(NLO) can be found in [30, 31]. In the quenched approx- 
imation of QCD, an effective low-energy chiral theory 
is formally obtained if an additional expansion in 1/N C , 
where N c is the number of colours, is carried out together 
with the usual one in quark masses and momenta [32, 33]. 
Here we adopt the pragmatic assumption that quenched 
ChPT describes the low-energy regime of quenched QCD 
in certain ranges of kincmatical scales at fixed N c . Corre- 
lation functions can be parametrized in terms of effective 
coupling constants, the latter being defined as the cou- 
plings that appear in the Lagrangian of the effective the- 
ory. For quark masses light enough to be in the e-regime 
of quenched QCD [34-36], the ratios corresponding to 
Eq. (5) in NLO ChPT, in a volume V = T x L 3 and at 
fixed topological charge v, are [7, 37] 

K±{x ,y ) = 1 ± ^ |/3i (|) 3/2 - /coo | , (10) 

where F is the pseudoscalar decay constant in the chi- 
ral limit, and the shape coefficients f3\ and fcoo can be 
found in Ref. [7]. Remarkably, the r.h.s. of Eq. (10) is 
determined once F is known, and it turns out to be in- 
dependent of v and the quark mass. When the quark 
masses are heavier and reach the so-called p-regime of 
QCD, the corresponding ratios are given by [31] 

M 2 / M 2 \ 
K (x ,y ) = 1 T 3 ( 47rF ) 2 log (jr) ± £0eo, yo) , (11) 

where M is the pseudoscalar meson mass at LO in ChPT, 
and fC(xo, yo) accounts for leading-order finite- volume ef- 
fects and can be found in Ref. [31]. The LECs gf can be 
extracted by requiring that 

gtK ± {x ,y ) = ktR ± (x ,yo) (12) 

for values of quark masses, volumes, x and yo, where 
quenched ChPT is expected to parametrize well the cor- 
relation functions. 

LATTICE COMPUTATION 

The numerical computation is performed by gener- 
ating gauge configurations with the Wilson action and 
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am 


aMp 


bare 


j^— , bare 


(R+ -i?-) baro 


e-regimc 


0.002 
0.003 




0.600(43) 
0.603(41) 


2.42(13) 
2.40(12) 


1.45(15) 
1.44(14) 


p-regime 


0.020 
0.030 
0.040 
0.060 


0.1960(28) 
0.2302(25) 
0.2598(24) 
0.3110(24) 


0.654(40) 
0.691(33) 
0.723(31) 
0.772(30) 


2.20(12) 
1.93(9) 
1.75(8) 
1.51(7) 


1.44(12) 
1.33(9) 
1.26(8) 
1.17(8) 



TABLE I: Results for aM P and ,R±' barc as obtained from 746 
and 197 gauge configurations in the e and p regimes, respec- 
tively. 



periodic boundary conditions by standard Monte Carlo 
techniques. The topological charge and the quark prop- 
agators are computed following Ref. [38]. The statisti- 
cal variance of the estimates of correlation functions has 
been reduced by implementing a generalization of the 
low-mode averaging technique proposed in [39], which 
turns out to be essential to get a signal for the lighter 
quark masses. The lattice has a bare coupling constant 
(3 = 6/.9q = 5.8485, which corresponds to a lattice spac- 
ing a ~ 0.12 fin, and a volume of Fa~ 4 = 16 3 x 32. The 
list of simulated bare quark masses, together with the 
corresponding results for pion masses and unrenormal- 
ized ratios R±< barc = Z 2 J R ± /Z±, are reported in Table I. 
Further technical details will be provided in a forthcom- 
ing publication. 

The values in Table I show that R±' barc exhibit a 
pronounced mass dependence, which is more marked in 
i?~ ,bare . We have explored several fit strategies, at- 
tempting to minimize the systematic uncertainties due 
to neglected higher orders in ChPT. The structure of 
Eqs. (10) and (11) indeed suggests that it is possible to 
cancel large NLO ChPT corrections by constructing suit- 
able combinations of R ± - barc . We observe that the prod- 
uct g^ is very robust with respect to the details of the 
fit strategy. The simplest way to extract this quantity is 
from a fit to the combination (i? + i?.~) baro , where NLO 
ChPT corrections cancel in the limit m —> 0. We obtain 



-\ba 



1.47(12) 



(13) 



To extract g+> bal ° anc ] g 1 ' hmc separately we then fit 
R+> baro to NLO ChPT, taking the value of F from a fit 
to the two-point functions as in Ref. [39] and the bare £ 
from Ref. [40]. Putting the result together with Eq. (13) 
we get 



9x 



0.63(4)(8), gl 



-,ba 



2.33(11)(30), (14) 



where the first error is statistical and the second is an es- 
timate of the systematic uncertainty from the spread of 
the central values obtained from fits to different quanti- 
ties and/or mass intervals. The physical LECs are given 
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FIG. 1: Mass dependence of R ± - h ^ c and {R+ ■ RT 
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(15) 



where kf 1 are the renormalization group-invariant (RGI) 
Wilson coefficients [1-4, 7]. The RGI quantities 



jj±,RGI _ ^,±,RGI 



1 w- 



--rlMl 



(16) 



at the pseudoscalar mass r^M^ = 1.5736 are taken from 
Refs. [41-44], and ro is a low-energy reference scale widely 
used in quenched QCD computations [45]. This proce- 
dure, analogous to the one proposed for the scalar den- 
sity in Ref. [46], provides values of the LECs that are 
non-perturbatively renormalized, as explained in detail 
in Ref. [44]. 

PHYSICS DISCUSSION 

By using the non-perturbative renormalization factors 
in Ref. [44] 



R+ 



RGI 



1.15(12) 



R 



.RGI 



i?". ba 



0.56(6), 



(17) 

and the perturbative values = 0.708 and k 1 = 1.978 
(see Ref. [7]), we obtain our final results 

g+ =0.51(9), 9l " = 2.6(5), g+g^ = 1.2(2) . (18) 

A solid estimate of discretization effects would require 
simulations at several lattice spacings, which is beyond 
the scope of this exploratory study. However, computa- 
tions of R^ at different lattice spacings and for masses 
close to m s /2 [7, 47] indicate that discretization effects 
may be smaller than the errors quoted above. It is also in- 
teresting to note that quenched computations of various 
physical quantities carried out with Neuberger fermions 
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show small discretization effects at the lattice spacing of 
our simulations [48, 49]. 

The values of gf in Eq. (18) are the main results of this 
paper. They reveal a clear hierarchy between the low- 
energy constants, ^ g^, which implies the presence of 
a Al =1/2 enhancement in the GIM- limit of (quenched) 
QCD. The strong mass dependence of R^' barc in Fig. 1 
indicates that an extrapolation of data around or above 
the physical kaon mass to the chiral limit is probably 
subject to large systematic uncertainties. 

When the charm mass m c is sufficiently heavier than 
the three light-quark masses, the chiral effective theory 
has a three-flavour SU(3) symmetry and the LO AS = 1 
effective Hamiltonian has two unknown LECs, 527 and 
g 8 - In our strategy these LECs are considered functions 
of the charm mass, and our normalizations are such that 2 

327(0) =g+, g s (0)=g- + ?t. (1 9 ) 



The values of gn^rfic) and g$(m c ) can be estimated at 
the physical value of the charm mass m c by matching 
the LO CHPT expressions with the experimental results 
for \Aq\ and l^l- The result is 



|<^ p (m c )|~0.50, 



\g? p (m c )\~10.5. 



(20) 



These estimates are, of course, affected by systematic er- 
rors due to higher-order ChPT contributions [51]. Keep- 
ing this in mind, the value of g^i^c) is in good agree- 
ment with our result. Since 327 is expected to have 
a mild dependence on the charm-quark mass (only via 
the fermion determinant in the effective gluonic action), 
and barring accidental cancellations among quenching ef- 
fects and higher-order ChPT corrections, this agreement 
points to the fact that higher-order ChPT corrections in 
\A 2 \ may be relatively small. Our value for gs(0) differs 
by roughly a factor of 4 from (?g Xp (m c ) given in Eq. (20). 
Apart from possible large quenching artefacts, our result 
suggests that the charm mass dependence and/or higher- 
order effects in ChPT are large for \A \. Indeed in this 
case penguin contractions, which are absent in the GIM 
limit, can be responsible for a large charm-mass depen- 
dence in g s , a dependence that can be studied in the next 
step of our strategy [7, 50]. 
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